Exciton and biexciton energies in bilayer systems 
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We report calculations of the energies of excitons and biexcitons in ideal two-dimensional bilayer 
systems within the effective-mass approximation with isotropic electron and hole masses. The exci- 
ton energies are obtained by a simple numerical integration technique, while the biexciton energies 
are obtained from diffusion quantum Monte Carlo calculations. The exciton binding energy decays 
as the inverse of the separation of the layers, while the binding energy of the biexciton with respect 
to dissociation into two separate excitons decays exponentially. 

PACS numbers: 02.70.Ss, 71.35.Cc, 78.67.De 
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Bound excitons and biexcitons have been observed in semiconductors under a variety of conditions. In this pa- 
per we consider bilayer systems in which the electrons are spatially separated from the holes, leading to what are 
known as "indirect excitons." Such systems have been realized in double-quantum-well structures under an applied 
perpendicular electric field, which serves to confine electrons in one well and holes in the otherjii^ The possibility of 
Bose-Einstein condensation of excitons in such structures has recently aroused much interestii^i^i^ and there is a need 
for a deeper understanding of the processes which laser excitation initiates in these systems. In this paper we consider 
an aspect of excitations in coupled quantum wells which may be relevant to experiments on coupled quantum wells - 
the energetics of biexcitons in bilayer systems. 

The effective-mass approximation with isotropic electron and hole masses gives a simple description of excitons and 
biexcitons which has been applied to many systems. This model is highly idealized, and effects due to anisotropic 
masses, non-parabolic bands, and finite well widths and depths will be significant. The model is, however, simple 
enough to be solved to very high accuracy, providing benchmark results while still permitting comparisons with 
experimental data. We have calculated exciton and biexciton energies within the effective-mass approximation for a 
system consisting of ideal two-dimensional electron and hole layers separated by a distance d. 

An exciton in an ideal two-dimensional bilayer geometry is described by the Schrodinger equation 
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where rue and are the electron and hole masses, respectively, and e is the static dielectric constant of the material. 
In the following, energies are given in terms of the exciton Rydberg, Ry* = /Xe;ie''/(2(47reoe)^ft^), and lengths are 
given in terms of the exciton Bohr radius, a*g — Aweoeh'^ / (/ie/ie^), where fieh = niemh/ {me + rah) is the reduced mass. 

Eq. ID) may be simplified by transforming into the center-of-mass frame and separating the variables in cylindrical 
polar coordinates. For the zero-angular-momcntum states we obtain 
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where r is the in-plane component of the electron-hole separation. Eq. ([2]) may be solved analytically when d = 0, 
which gives a ground state wave function of <I>(r) oc exp[— 2r] and an energy of Ex — —4 Ry* . For d > we solved 
Eq. ([2]) using a standard Runge-Kutta numerical integration technique. The exciton energy is plotted as a function 
of d in Fig. [TJ The energy takes its minimum value of —4 Ry* at d = 0, while at small separations the energy varies 
linearly with d and at large separations it varies as \jd. The results shown in Fig. [1] may be fitted to the expression 
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where A = 154.363, B = 648.9, C = 225.005, D = 46.4263, E = 384.976, F = 628.158, and G = 129.672. 
expression gives a maximum error of less than 0.0028 Ry* in the range < d < 10 a|j. 
The Schrodinger equation for the biexciton is 
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FIG. 1: Exciton energy as a function of the separation of the electron and hole layers. 

where 1 and 2 denote the electron coordinates, a and b denote the hole coordinates, ri2 = |ri — r2|, ria = 
\/\ri — Top + (P^ etc., and a = rrie/mh- When expressed in units of Ry* , Ex is a function only of d, but the 
biexciton energy, Exx, is a function of both d and a. Eq. ([4]) does not separate in cylindrical polar coordinates and 
we have to solve the many-body problem. For this purpose we have used the diffusion quantum Monte Carlo (DMC) 
method, which is a stochastic projector technique for solving the imaginary-time many-body Schrodinger equation.^ 
In the ground state of the biexciton the electrons have opposite spins and the holes have opposite spins, so the spatial 
part of the wave function is node-less. The DMC method is exact in principle for node-less wave functions, and 
although there are biases due to the use of finite time steps and populations of walkers, these can be made negligible 
for small systems such as this. 

The sampling within DMC is guided by an approximate wave function which must be sufficiently accurate to give 
low statistical noise and to keep the biases small. The form of our approximate wave function was guided by the 
symmetries of the problem and the long- and short-distance behavior. The system composed of two separated bound 
excitons is always more stable than one consisting of four unbound charges. Therefore we expect the wave function 
to be exponentially small when all four particles are far apart. When one of the particles is far from the other three 
we expect the wave function to be exponentially small because the single charge will be attracted to the other three. 
Likewise we expect that the part of the wave function corresponding to one bound exciton and a free electron and hole 
is exponentially small. When d is large we expect the system to consist essentially of two separated excitons, and the 
form of the approximate wave function must allow for this possibility. The short range behavior of the wave function 
is fixed by the Kato cusp conditions^, which ensure that the divergences in the potential and kinetic energies cancel 
when two particles are coincident. The biexciton wave function, '5, should be unaltered by exchange of (i) the two 
electron coordinates, or (ii) the two hole coordinates, i.e., ^'(ri, r2, Tq, r^) = ^'(r2, ri, Tq, rf,) = ^'(ri, r2, r^, Tq), and 
when the electron and hole masses are equal, should have the additional electron-hole symmetry ^'(ri, r2, r^, r^) = 
*{ra,r6,ri,r2). 

The binding between the excitons is expected to be small compared with the binding within an exciton. We 
therefore write the wave function as an appropriately-symmetrized product of two exciton wave functions, which is 
then multiplied by a Jastrow function containing electron-electron and hole-hole terms. We use the following form. 
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which satisfies all of the above conditions, 
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where Ci-Cg are parameters. This form is similar to that used by Lee et al7 for the three-dimensional biexciton, 
although our form has more flexibility in the electron-hole part. We require C2, C4, cy, cg > so that ^ is well-behaved, 
and Cg, C8 < so that eh decays when the electrons and holes are far apart. Eq. ([5]) describes two separated excitons 
when either cg or cg go to zero. 

The values of the parameters ci and C3 were fixed by the electron-electron and hole-hole Kato cusp conditions.— The 
value of C5 was fixed by the electron-hole cusp condition when d = 0, while for d > there should be no electron-hole 
cusp and so we set C5 = 0. When a = \ electron-hole symmetry requires that ci = C3 and C2 — c/^. The optimal 
values of the remaining variable parameters were obtained by minimizing the variance of the variational energy ..^'^^ 

We calculated the energy of the biexciton, Exx, as a function of d for cr = 1 and 2. Tests indicated that the timestep 
and population control errors in the DMC results were negligible. All the QMC calculations were performed using the 
CASINO code.^" The biexciton binding energy with respect to dissociation into two separate excitons, Eb = 2Ex — Exx, 
is plotted as a function of d in Fig. [5] for a — 1 and 2. For small d and cr = 1, i?f, is close to the value obtained in 
earlier calculations for d = of 0.771 Ry*!^ Et goes to zero at large d much more rapidly than Ex itself because the 
electron-electron and hole-hole repulsions dominate the electron-hole attraction in the biexciton at large d, tending 
to unbind the biexciton. Examination of the biexciton wave function shows that two separated excitons are formed 
at large d. The behavior of Eb at large d is reasonably well represented by a simple exponential form, 



Eb = a exp[— /3d] , 
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where a = 0.67573 and f3 = 15.023 for = 1, and a = 0.71714 and (3 = 14.187 for ct = 2. 




FIG. 2: The binding energy of the biexciton, Eb = 2Ex — Exx, as a function of the separation of the electron and hole layers, 
d, for (7=1 and 2. The error bars are smaller than the thickness of the lines. 



As a simple example of the use of these results we estimate the biexciton binding energy, Eb, in the experiments 
of Butov et alX, who studied a system of two 80 A-wide GaAs quantum wells separated by a 40 A-wide barrier of 
Alo.33Gao.67As. The electron mass in GaAs is rrif, = 0.067mQ, while the heavy-hole mass should be reduced from its 
bulk value of 0.45mo by confinement effects, and for simplicity we take a value of rrih = 0.134mo. This gives a 2:1 
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mass ratio, although the results are not sensitive to the precise value of m^. Using a dielectric constant appropriate 
to GaAs of 13.2 we find = 156 A and Ry* = 3.5 meV. Ef, is sensitive to the value of d, and therefore we fix its 
value such that we reproduce the large-field exciton binding energy of 4 meV (1.14 -Ry*) calculated for this structure 
by Szymanska and Littlewood^^, who used a realistic description of the finite well widths and depths. Using Eq. ^ 
we find d — 0.64 a*g (100 A), which is a very reasonable value as it lies between the experimental barrier width of 
40 A and the distance between the centers of the wells of 120 Aji Substituting d — 0.64 a*Q into Eq. ^ and using 
the parameters for cr = 2 we obtain Ei, = 8.2 x 10~^ Ry* (2.9 x 10"''' meV). This model predicts an extremely small 
biexciton binding energy, which is unlikely to lead to measurable effects. 

In conclusion, we have calculated the energies of excitons and biexcitons in ideal two-dimensional bilayer systems 
within the effective-mass approximation with isotropic electron and hole masses. The exciton binding energy decays 
as the inverse of the layer separation, while the biexciton binding energy with respect to dissociation into two separate 
excitons decays exponentially. This model predicts that the biexciton binding energy in the experiments of Butov et 
alX is extremely small. 

We thank Peter Littlewood for useful discussions. NDD and RJN acknowledge financial support from the Engi- 
neering and Physical Sciences Research Council (EPSRC), UK. 
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